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We present a theoretical study of the optical conductivity of graphene with magnetic and non¬ 
magnetic adatoms. First, by introducing alternating potential in a pure graphene, we demonstrate 
a gap formation in the density of states and the corresponding optical conductivity. We highlight 
the distinction between such a gap formation and the so-called Pauli blocking effect. Next, we apply 
this idea to graphene with adatoms by introducing magnetic interactions between the carrier spins 
and the spins of the adatoms. Exploring various possible ground-state spin configurations of the 
adatoms, we find that antiferromagnetic configuration yields the lowest total electronic energy, and 
is the only configuration that forms a gap. Furthermore, we analyze four different circumstances 
leading to similar gaplike structures and propose a means to interpret the magneticity and the pos¬ 
sible orderings of the adatoms on graphene solely from the optical conductivity data. We apply this 
analysis to the recently reported experimental data of oxygenated graphene. 

PACS numbers: 73.22.Pr, 78.67.Wj 


I. INTRODUCTION 


II. REAL GAP VS PAULI BLOCKING 


Setting adatoms on graphene has been a promising 
means to tailor electronic properties of graphene, such 
as transport, optical, and magnetic properties^P— . In 
particular, the idea of generating a controllable band 
gap in graphene while at the same time adding mag¬ 
netic moments through the adatoms has spawned a broad 
interest to make graphene a magnetic semiconductor. 
Theoretically, it has been predicted that nonmagnetic 
atoms could become magnetic when put as adatoms on 
graphene^. An experimental evidence of this is found in 
fluorinated graphene^. However, to make the magnetic 
moments align to form a ferromagnetic order is still a 
great challenge to both experimentalists and theorists^. 
Meanwhile, available theories based on Ruderman-Kittel- 
Kasuya-Yosida (RKKY) mechanism predict that a pair 
of neighboring magnetic moments of the adatoms on 
graphene prefer to align antiferromagnetically^iS. Indi¬ 
cation of such antiferromagnetic alignments has been re¬ 
ported recently in fluorinated graphene, for which spin- 
half paramagnetism is observed with only a very small 
portion of magnetic moments of the fluorine adatoms 
contribute to the magnetization^. 

In this paper, we aim at showing the intimate connec¬ 
tion between magneticity of the adatoms on graphene, 
their magnetic orders, i.e. antiferroinagnetic (AFM), fer- 
omagnetic (FM), and paramagnetic (PM), and their cor¬ 
responding optical responses. In particular, we address 
how a band gap can form in the density of states (DOS) 
of this system, and how such a gap reveals in the op¬ 
tical conductivity ((Ti(w)). We analyze several circum¬ 
stances leading to similar gaplike structures and compare 
our results with the recently reported experimental data 
of (Ti(a;) of graphene with oxygen adatomsiS at various 
adatom concentrations. 


Before we go through the more realistic model for 
graphene with magnetic adatoms, let us consider a toy 
model of graphene with some artificial potensial V(—V) 
at each sublattice A{B). Such an alternating potential 
has been used in the discussion on excitonic mass gen¬ 
eration in grapheneii. Using the tight-binding approx¬ 
imation defined on a hexagonal lattice as shown as the 
honeycomb net in Fig. [TJ the Hamiltonian for graphene 
with alternating potential (neglecting electron-electron 
interactions) can be written in k-space as 
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In Eq. ([T]) aj^^(ak,o-) creates(annihilates) an electron at 
sublattice A, while 6]^ ^.(fck.o-) creates (annihilates) an elec¬ 
tron at sublattice B, with momentum k and spin a. :/3(k) 
is the nearest-neighbor hopping energy dispersion, with t 
being the nearest neighbor hopping integral and a being 
the lattice constant. 

Since the Hamiltonian contains no spin-dependent in¬ 
teractions, we can define the corresponding Green’s func¬ 
tion matrix as 


[G(k,2)] 
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with z = w ± iO’*', and oj real frequency variable. Here, 
we set h = 1 such that frequency and energy have the 





FIG. 1: (Color online) Illustration of graphene in the presence 
of adatoms. The honeycomb net represents a graphene lat¬ 
tice. The black filled and opaque circles represent adatoms sit¬ 
ting at sublattice A and B, respectively. Parallelograms with 
different colors and line styles represent supercells with dif¬ 
ferent configurations of adatom occupation: black-solid lines 
(dotted-red lines) for a supercell containing only an adatom at 
sublattice A{B), dashed-green lines for a supercell containing 
adatoms at both sublattices, and dot-dashed-blue lines for a 
supercell with no adatoms. 



FIG. 2: (Color online) Comparison between optical conduc¬ 
tivity of (a) bare graphene, (b) graphene with an alternat¬ 
ing potential, and (c) charged graphene. The magnitude 
of the alternating potential and the chemical potential are 
taken to be as follows: V = fj, = 0 for the bare graphene, 
V = 0.5 eV, /i = 0 for graphene with alternating potential, 
and V — 0, fi = 0.5 eV for charged graphene, respectively. 
The insets show the DOS for the three cases where the small 
vertical arrowed line in each inset indicates the position of the 
chemical potential. 


same unit. We can then use this Green’s function matrix 
to compute the DOS 

+ ( 4 ) 
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and the optical conductivity tensor derived from the 
Kubo formula, as previously used in Refs. [mi, 
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with [fA(k)] = d[H{'k)]/dk\ as the Cartesian component 
of the velocity matrix, 

[Aik, 0] = ^{[G(k, i/ + t0+)] - [G(k, r. - *0+)]}, (6) 

as the spectral function matrix, /(r^,/r,r) as the Fermi 
distribution function, u and v as the real frequency vari¬ 
ables, n as the chemical potential, and T as the tempera¬ 
ture. Here, the transverse component Uxyito) = ay^iuj) ss 
0, hence we shall address only the longitudinal compo¬ 
nent (Ti(w) = ffxxi^) ~ o’yy(w)- Furthermore, we are 
interested in the ground-state properties, i.e., at T = 0. 

Let us compare the results for bare graphene {V = 

= 0), graphene with alternating potential {V = 0.5 
eV, ^ = 0), and charged graphene (H = 0, /x = 0.5 eV), 
that are shown in Fig. [2j cri(a;) for the bare graphene 
[see the black curve (a) in the main panel] has a hnite 
value at the dc limit which is the universal value tro = 
7re^/2/i and has a peak associated with the transition 


between the two Van Hove singlarity points of the DOS. 
Here, we take t = 2.35 eV such that the position of the 
peak agrees with the experimental data of 171 ( 0 ;) we will 
address later. As we introduce the alternating potential, 
the gap of size A r; 2V = I eV forms in the DOS [see 
inset (b)]. Corresponding to the gap formation in the 
DOS, cri(a;) reveals a complete loss of spectral weight 
from o = 0 to O' « A [see red curve (b) in the main panel]. 
To compensate the /-sum rule, this loss of spectral weight 
is accompanied with the appearance of a ’’hump” at the 
gap edge. 

Now, let us see the situation for the charged graphene. 
Since the electrons are assummed noninteracting, shift¬ 
ing the chemical potential (/x) does not alter the profile 
of DOS compared to that of the bare graphene. How¬ 
ever, the shift of the position of fi significantly affects 
the profile of aiiur). Here we see a similarity between 
the prohle of (Ti(a;) for graphene with a gap and charged 
graphene without a gap [compare curves (b) and (c) of 
cri(a;)] in the energy region 0 < w < I eV but not quite 
close to 0. In this region the green curve (c) also shows 
a vanishing spectral weight akin to the formation of gap. 
This is known to be due to the Pauli blocking effecti^. 
It is important, however, to note the difference. Charged 
graphene has a finite DOS at /x, hence it is a metal rather 
than an insulator. This metallic characteristics is re¬ 
vealed by the formation of a sharp Drude peak at w = 0 
(dc conductivity). Therefore, an optical probe of 
not close enough to the dc limit may not be able to distin¬ 
guish whether the system is gapped or only manifesting 
the Pauli blocking effect. 

Note that if /x is too small, away from the Drude 
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response cri(a;) may not completely vanish before it 
rises back. Such a situation reveals in graphene under 
a gate voltage as high as 71 volta^^, from which we 
may infer that 2fi Ri 5000 cm“^, ot ^ ~ 0.3 eV. This 
situation may be explained as the following. In a 
Pauli-blocked system of graphene, the Drude peak has 
a width, and hence a tail, which samples the finite 
values of the density of states closely below and above 
/i. If fi is not big enough this tail effect would produce 
some residual conductivity in the energy region between 
0 and 2/r. Theoretical studies on graphene with a 
shifted chemical potential that incorporate interactions 
among electrons and between electrons and other quasi¬ 
particles suggest additional broadening of the Drude 
peak and renormalization of the gap edge associated 
with the interband transitio n^^d'^ . These factors, in 
turn, enhance the residual conductivity mentioned above. 


In the second line of Eq. 0 ) denotes the spin of 

an adatom on top of site i{j) of sublattice A{B), whereas 
SAi(s_Bj) denotes the spin of an electron correspondingly. 
Jh is Hund’s coupling between an electron spin and spin 
of an adatom at the same site. It is customary to treat 
the local spins in the above Hamiltonian classically. Such 
a problem would then be easily solved within the dy¬ 
namical mean-field theory (DMFT) ^°d^ as similarly per¬ 
formed in Ref. 

To start the DMFT algorithm we define the Green’s 
function matrix through the Dyson equation, 

[G(k,z)]-i = [Go(k,z)]-i-[E(z)], (8) 

with 


III. MODEL FOR GRAPHENE WITH 
ADATOMS 


In what follows we discuss a model for graphene with 
magnetic adatoms. Despite the fact that adatoms on 
graphene may tend to form clustersi^ or ordered statesi^, 
for simplicity we assume that they are distributed ran¬ 
domly over the graphene sheet and each adatom rests 
on top of a carbon atom. Suppose we have x portion 
of adatoms out of N lattice sites. Thus, statistically we 
expect to have a:(I — x) portion of supercells containing 
only an adatom at sublattice A, the x(I —x) portion con¬ 
taining only an adatom at sublattice B, the x^ portion 
containing adatoms at both sublattices A and B, and the 
(I — x)^ portion with no adatoms [see the illustration in 
Fig. d]. 

Considering that each adatom carries a magnetic mo¬ 
ment corresponding to its spin, we introduce local mag¬ 
netic interactions between the carrier spins and the spins 
of the adatoms such that the full Hamiltonian reads 
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the unperturbed Green’s function matrix, and [S(z)] as 
some initial guess of self-energy matrix. Then, we av¬ 
erage [G(k,z)] over the first Brillouin zone of graphene, 
[G{z)\ = ^k[^(kj •2^)]- The true [E(z)] is to be solved 
through a self-consistency process. 

The mean-field Green’s function matrix is extracted 
through [g{z)] = [[G(z)]-i -f [E(z)]]-b This is then 
used to compute the local interacting Green’s function 
matrix for a supercell with: an adatom at sublattice A 
only, [Ga[z^6at4>a)\\ an adatom at sublattice B only, 
[Gb(z, 0s, 0 b)]; and adatoms at both sublattices A and 
B, [Gab{z,0a,4'a,(^b,4>b)]- Each inverse of these local 
interacting Green’s function matrices is obtained by sub- 
stracting its corresponding local interaction matrix, 
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from [Q{z)]~^. The three local interacting Green’s function matrices, corresponding to the supercells with different 
adatom occupations are to be averaged over all possible directions of adatom spins Sa and Sb in their corresponding 
solid angles Ha and Hb, 


[Ga{z)] = J dHAPAi9A, 4’A)[GAiz, 9 a, 4>a)], 

[Gb(z)] = y C?nBPB(0B, (^b)[Gb(2:, 0B, <^b)], (10) 

[Gab(2:)] = J dHAdHBPABi9A,4‘A,9B,4’B)[GAB{z,9A,4’A,9B,4>B)]- 

Here, Pa{9a,4>a) and Pb{9b,4>b) are assumed constant for all angles since Sa or Sb standing alone in a supercell 
is assumed to be free to point to any direction. Whereas, Pab{9a, 4>a, 9b, 4>b) depends on the assumed ground-state 
magnetic configuration of the adatoms, namely 


6{9a)S{9b - 

Pab{9a,4>a,9b,4>b) oc < S{9a)S{9b) 

constant 


The average interacting Green’s function matrix is 
then obtained by further averaging over supercells with 
different adatom occupations, 

[G{z)]ave = - a;)[GA(-2)] -I- x{l - x)[Gb(z)] 

+x^[Gab{z)] +{1 - x^lGiz)]. (11) 

Finally, the corrected self-energy matrix is extracted 
from the Dyson equation [S(z)] = [G{z)]~^ — [G{z)]~^^, 
which is then used to start the new iteration by feeding 
it back into Eq. (|S]). This self-consistency process is con¬ 
tinued until [^(z)] is converged. Once the convergency is 
achieved, we compute the DOS and fTi(w) according to 
the definions given by Eq. (U) and Eq. ®, respectively, 
except that now the summation over the spin index 
is not needed as the spin dependence has been taken 
explicitly into the definition of the matrices. 

IV. DISCUSSION ON VARIOUS MAGNETIC 
GROUND STATES 

Eigure [3] displays our calculation results of optical 
conductivities of graphenes with 100% adatom cover¬ 
age (x = 1) for various possible magnetic configurations. 
Note that, despite the fact that the calculations we per¬ 
form are for the ground states, i.e., T = 0, the results 
should still be relevant to compare with finite temper¬ 
ature experimental data as long as the temperature is 
well below any magnetic ordering transition tempera¬ 
ture. For these calculations we choose the parameter 


tt) for the AFM configuration, 
for the FM configuration, 
for the PM configuration. 

r 



FIG. 3: (Color online) Comparison of optical conductivities 
of graphenes with 100% adatom coverage for various possible 
magnetic configurations. Main panel shows the optical con- 
ductiviy for (a) antiferromagnetic (AFM), (b) ferromagnetic 
(FM), and (c) paramagnetic (PM) ground-states, respectively. 
The insets show their corresponding DOS profiles. 


values Sa = Sb = 1/2, Jh = 1-2 eV, and t = 2.35 eV. 
The main panel shows the distinction among (Ti(a;)’s for 
AFM, FM, and PM configurations. From the DOS pro¬ 
files shown in the insets, it is clearly seen that a gap forms 
only for the AFM ground state. We also compute the to¬ 
tal electronic energy, u!p(u!)duj for each of the three 
possible ground-state configurations and confirm that the 

















AFM configuration gives the lowest energy. The physical 
concept of the formation of the gap in the AFM configu¬ 
ration is essentially that the AFM configuration induces 
alternating potentials as discussed previously above. 

At this point, we would like to address the experimen¬ 
tal data of fluorinated graphene presented in Ref. @ 
in connection with our results. Because the vanishing 
contribution of the most magnetic moments is due 
to clustering of the fluorine adatoms, we argue that 
every pair of neighboring magnetic moments inside the 
cluster align antiferromagnetically, whereas only the 
magnetic moments at the edge have the freedom to 
change their orientations. To verify our prediction, an 
optical conductivity measurement is necessary for which 
we predict a profile of cri(a;), such as the black curve in 
Fig. [3] would result. 


V. ANALYSIS OF OPTICAL 
CONDUCTIVITY WITH SIMILAR GAPLIKE 
STRUCTURES 

In the following we analyze for four different 

circumstances leading to similar gaplike structures. Our 
aim is to propose a means to predict the magneticity and 
the possible orderings of the adatoms on graphene solely 
from the experimental optical conductivity data. As an 
example, we apply this analysis to the recently reported 
experimental data of oxygenated graphenei^ [replotted 
in panel (a) of Fig. |4]. Since, unfortunately, the experi¬ 
mental data of cri(a;) for photon energy below 0.55 eV are 
not yet available, we focus on analyzing the behavior of 
ai(u}) within the available energy range 0.55 — 5.3 eV). 
The result of this analysis may later be verified when the 
low energy data have become available. Note that our 
assumption on the possibility that oxygen adatoms pos¬ 
sess magnetic moments is based on several experimental 
hints2^“— suggesting that oxygen can be magnetic in cer¬ 
tain environments. 

The experimental data of cti(u;) displayed in Fig. |l](a) 
show that a gaplike structure gradually forms, whereas 
the 4.7-eV peak monotonically decreases its intensity as 
the system adsorbs more oxygen adatoms. To compare 
with our calculations, first, we assume that the adatoms 
are magnetic and form an AFM order and explore the 
system for various adatom concentrations (cc) as shown 
in panel (b). As x is increased, a real gap emerges ac¬ 
companied by the appearance of a hump at the gap edge. 
The 4.7-eV peak intensity decreases as x is increased up 
to 0.6 and increases back accompanied with a blue shift 
for a further increase in x. Apparently, the appearance 
of a hump and the nonmonotonic change in intensity of 
the 4.7-eV peak do not agree with those observed in the 
experimental data. 

Second, as shown in panel (c), we explore the system 
with FM Oder for various cc’s as performed for the AFM 
case. As x is increased, a gaplike structure with a Drude 
peak forms. Apart from the Drude peak, this behavior 



FIG. 4: (Color online) Comparison of between exper¬ 

imental data and calculation results, (a) The experimental 
data of oxygenated graphene for various times of oxygen ex¬ 
posure (replotted from Ref. [3). (b) and (c) Graphene sys¬ 
tems with FM- and AFM- ordered adatoms, respectively, for 
various adatom concentrations (x). (d) Graphene with non¬ 
magnetic adatoms for various chemical potential values (p). 
(e) Graphene with magnetic adatoms of the PM configuration 
with p increasing proportionally to x. The vertical red dashed 
line is to mark the energy value below which the experimental 
data are not available for comparison. 
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looks similar to that in the experimental data. The 4.7- 
eV peak also changes its intensity as x is increased up to 
0.6, but it increases back as x is increased further. Again, 
this nonmonotonic change in intensity of the 4.7-eV peak 
does not agree with that observed in the experimental 
data. 

Third, we consider nonmagnetic adatoms and assume 
that they supply some extra charge to graphene, such 
that the chemical potential (/r) is shifted proportionally 
to X. As shown in panel (d), shifting /r away from 0 leads 
to the formation of a gaplike structure with a Drude peak 
in the low-energy region. As far as for the energy region 
above 0.55 eV, this behavior looks similar to that in the 
experimental data. However, we notice that the 4.7-eV 
peak does not change as jx changes, contradictory to the 
experimental data. 

Last, we consider graphene with magnetic adatoms in 
PM configuration [see panel (e)]. Here we also assume 
that /r varies proportionally to x. By increasing x, hence 
also /r, we generate a gaplike structure with a Drude 
peak. Apart from the Drude peak, remarkably we see 
that the 4.7-eV peak monotonically decreases as x in¬ 
creases, consistent with that observed in the experimen¬ 
tal data. Thus, within the energy range of 0.55-5.3 eV, 
aiiio) of the oxygenated graphene seems to resemble that 
of graphene with the adatoms of the PM configuration. 

Now let us elaborate further our results in general. 
We have seen that a hump structure appears in (Ti(w) 
whenever the system forms a gap, suggesting that the 
hump may be a signature of the presence of a gap in 
a graphene sample. On the other hand, the absence of 
a hump of a gaplike structure in cri(a;) may be the sig¬ 
nature of the presence of a Drude peak or the Pauli- 
blocking effect. Meanwhile, by inspecting whether or not 
the 4.7-eV peak changes its intensity as a function of 
adatom concentration, we may predict that the adatoms 
posses magnetic moments if the peak intensity changes 
with variation of adatom concentration, otherwise, the 
adatoms are nonmagnetic. Interestingly, how the 4.7-eV 
peak changes its intensity as a function of adatom con¬ 
centration may signify a possible magnetic ordering of 
the adatoms. Namely, the adatoms may form: AFM or¬ 
der if the 4.7-eV peak shows a nonmonotonic change in 
its intensity, accompanied by a blue shift; FM order if the 
4.7-eV peak shows a nonmonotonic change in its inten¬ 
sity, without a blue shift; or PM configuration if the 4.7- 
eV peak shows a monotonic dicrease in its intensity; as 
the adatom concentration is increased. Such connections 
between the optical conductivity data and the magnetic- 
ity and the ordering of the adatoms, if true, may provide 
a promising means to deduce the magnetic characteris¬ 
tics of adatoms on graphene in case magnetic probes are 
not available. We therefore suggest that it be crucial to 
verify these connections by experiments with magnetic 
probes and/or optical conductivity down to very low en¬ 


ergy region close to the dc limit. 

Lastly, we comment on the little dissimilarity of our 
results as to compare with the experimental data of 
oxygenated graphene. Namely, the 4.7 eV peak appears 
to be much broader in the experimental data than that 
obtained in our calculations. We conjecture that this is 
due to some factors not taken into account in our model, 
possibly the effects of a substrate, or interplay between 
electron-electron Coulomb interactions and interactions 
between electrons spins and the spins of the adatoms. 
These demand further theoretical studies. 


V. CONCLUSION 

We have presented our theoretical study of graphene 
with magnetic and nonmagnetic adatoms. We start 
with an introductory study of systems of graphene 
with alternating potentials and a bare graphene with a 
shifted chemical potential. We demonstrate the optical 
conductivities corresponding to these systems. Both 
systems posses similar gaplike structures in cri(a;), but 
they differ in that only graphene with the alternating 
potential reveals a real gap, whereas graphene with 
the shifted chemical potential reveals a sharp Drude 
peak. Furthermore, we study a model for graphene with 
adatoms with magnetic interactions between electron 
spins and spins of the adatoms. We demonstrate the 
different profiles of cri(w) corresponding to AFM, FM, 
and PM adatom configurations. We show that only 
graphene with AFM adatom configurations possesses a 
gap, and this configuration gives the lowest electronic 
energy. Furthermore, we analyze four different circum¬ 
stances leading to similar gaplike structures in fTi(w) and 
propose connections between characteristics revealed in 
cri(a;) and the magneticity and the spin ordering of the 
adatoms. We compare our results with the experimental 
data of oxygenated graphene from Ref. [l^ . from which 
we conjecture that the adatoms are magnetic with the 
PM configuration. Our results show a path to modify 
optical properties as well as band gap of graphene via 
magnetic adatoms. 
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